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Abstract: The muon anomalous magnetic moment is one of the most precisely 
measured quantities in particle physics. Recent high precision measurements 
(0.54ppm) at Brookhavcn reveal a "discrepancy" by 3.2 standard deviations 
from the electroweak Standard Model which could be a hint for an unknown 
contribution from physics beyond the Standard Model. This triggered nu- 
merous speculations about the possible origin of the "missing piece". The 
remarkable 14-fold improvement of the previous CERN experiment, actually 
animated a multitude of new theoretical efforts which lead to a substantial 
improvement of the prediction of a M . The dominating uncertainty of the pre- 
diction, caused by strong interaction effects, could be reduced substantially, 
due to new hadronic cross section measurements in electron-positron annihi- 
lation at low energies. After an introduction and a brief description of the 
principle of the experiment, I present a major update and review the sta- 
tus of the theoretical prediction and discuss the role of the hadronic vacuum 
polarization effects and the hadronic light-by-light scattering contribution. 
Prospects for the future will be briefly discussed. As, in electroweak precision 
physics, the muon g — 2 shows the largest established deviation between the- 
ory and experiment at present, it will remain one of the hot topics for further 
investigations. 

1 Lepton magnetic moments 

The subject of our interest is the motion of a lepton in an external electromag- 
netic field under consideration of the full relativistic quantum behavior. The 
latter is controlled by the equations of motion of Quantum Electrodynamics 
(QED), which describes the interaction of charged leptons [l = e,fj,,r) with 
the photon (7) as an Abelian f7(l) e m gauge theory. QED is a quantum field 
theory (QFT) which emerges as a synthesis of quantum mechanics with special 
relativity. In our case an external electromagnetic field is added, specifically a 
constant homogeneous magnetic field B. For slowly varying fields the motion 
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is essentially determined by the generalized Pauli equation, which also serves 
as a basis for understanding the role of the magnetic moment of a lepton on 
the classical level. As we will see below, in the absence of electrical fields E the 
quantum correction miraculously may be subsumed in a single number, the 
anomalous magnetic moment a^, which is the result of relativistic quantum 
fluctuations, usually simply called radiative corrections. 

Charged leptons in first place interact with photons, and photonic radiative 
corrections can be calculated in QED, the interaction Lagrangian density of 
which is given by (e is the magnitude of the electron's charge) 

C^ D (x) = ej^(x) A,(x) , %Jx) = - Mx^Mx) , U) 

i 

where j£ m (x) is the electromagnetic current, ipe(x) the Dirac field describing 
the lepton £, 7 M the Dirac matrices and with a photon field A fl (x) exhibiting 
an external classical component and hence A^ — > A^ + A° xt . We are thus 
dealing with QED exhibiting an additional external field insertion "vertex" . 

Besides charge, spin, mass and lifetime, leptons have other very interesting 
static (classical) electromagnetic and weak properties like the magnetic and 
electric dipole moments. Classically the dipole moments can arise from either 
electrical charges or currents. A well known example is the circulating current, 
due to an orbiting particle with electric charge e and mass m, which exhibits 
a magnetic dipole moment fiL — r x v given by 

»l = ^L (2) 

Zmc 

where L = m r x v is the orbital angular momentum (r position, v velocity) . 
As we know, most elementary particles have intrinsic angular momentum, 
called spin, and in particular leptons like the electron arc Dirac fcrmions of 
spin i . Spin is directly responsible for the intrinsic magnetic moment of any 
spinning particle. The fundamental relation which defines the "^-factor" or 
the magnetic moment is 

eh 

fi = gi S , S the spin vector. (3) 

2mgc 

For leptons, the Dirac theory predicts gi = 2 [I], unexpectedly, twice the 
value 5 = 1 known to be associated with orbital angular momentum. It took 
about 20 years of experimental efforts to establish that the electrons mag- 
netic moment actually exceeds 2 by about 0.12%, the first clear indication of 
the existence of an "anomalous" contribution to the magnetic moment [2] . In 
general, the anomalous magnetic moment of a lepton is related to the gyro- 
magnetic ratio by 

at = fit/ fi B - 1 = , (£ = e,n,r) (4) 
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where hb is the Bohr magneton which has the value 

eh 

[is = -z = 5.788381804(39) x 1CT 11 MeVT -1 . (5) 

2m e c 

Formally, the anomalous magnetic moment is given by a form factor, de- 
fined by the matrix element 

(*-(pO|j£n(o)|r(p)> 

where \£~(p)) is a lepton state of momentum p. The relativistically covariant 
decomposition of the matrix element reads 



(-ie)u(p>) T»F E tf)+iZ£FF M tf) u(p) 



with q = p' — p and where u(p) denotes a Dirac spinor, the relativistic wave 
function of a free lepton, a classical solution of the Dirac equation (7^^ — 
to) u(p) — 0. Fft{q 2 ) is the electric charge or Dirac form factor and Fm(? 2 ) 
is the magnetic or Pauli form factor. Note that the matrix a^ v = ^[7 A ',7 I/ ] 
represents the spin 1/2 angular momentum tensor. In the static (classical) 
limit q 2 — > we have 

F E (0)=1 ; F M (0) = a„ (6) 

where the first relation is the charge normalization condition, which must be 
satisfied by the electrical form factor, while the second relation defines the 
anomalous magnetic moment. a M is a finite prediction in any renormalizablc 
QFT: QED, the Standard Model (SM) or any renormalizable extension of it. 

By end of the 1940's the breakthrough in understanding and handling 
renormalization of QED (Tomonaga, Schwinger, Feynman, and others) had 
made unambiguous predictions of higher order effects possible, and in particu- 
lar of the leading (one-loop diagram) contribution to the anomalous magnetic 
moment 

af D W = ^,(£ = e,M,r) (7) 

by Schwinger in 1948 [3]. This contribution is due to quantum fluctuations 
via virtual photon- lepton interactions and in QED is universal for all lcptons. 
At higher orders, in the perturbative expansion^, other effects come into play: 
strong interaction, weak interaction, both included in the SM, as well as yet 
unknown physics which would contribute to the anomalous magnetic moment. 

In fact, shortly before Schwinger's QED prediction, Kusch and Foley in 
1948 established the existence of the electron "anomaly" g e — 2 (1.00119 ± 
0.00005), a 1.2 per mill deviation from the value 2 predicted by Dirac in 1928. 



1 which is equivalent to the loop-expansion, referring to the number of closed loops 
in corresponding Feynman diagrams. 
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We now turn to the muon. A muon looks like a copy of an electron, which 
at first sight is just much heavier ra fi /m e ~ 200, however, unlike the electron 
it is unstable and its lifetime is actually rather short. The decay proceeds by 
weak charged current interaction into an electron and two neutrinos. 

The muon is very interesting for the following reason: quantum fluctua- 
tions due to heavier particles or contributions from higher energy scales are 
proportional to 

where M may be 

- the mass of a heavier SM particle, or 

- the mass of a hypothetical heavy state beyond the SM, or 

- an energy scale or an ultraviolet cut-off where the SM ceases to be valid. 

On the one hand, this means that the heavier the new state or scale the harder 
it is to see (it decouples as M — ► oo). Typically the best sensitivity we have for 
nearby new physics, which has not yet been discovered by other experiments. 
On the other hand, the sensitivity to "new physics" grows quadratically with 
the mass of the lepton, which means that the interesting effects are magnified 
in a M relative to a e by a factor (m^/m e ) 2 ~ 4 x 10 4 . This is what makes 
the anomalous magnetic moment of the muon the predestinated "monitor for 
new physics" or, if no deviation is found it may provide severe constraints to 
physics beyond the 

In contrast, a e is relatively insensitive to unknown physics and can be 
predicted very precisely, and therefore it presently provides the most precise 
determination of the fine structure constant a = e 2 /4ir. 

What makes the muon so special for what concerns its anomalous magnetic 
moment? 

• Most interesting is the enhanced high sensitivity of a M to all kind of inter- 
esting physics effects. 

• Both experimentally and theoretically a M is a "clean" observable, i.e., it 
can be measured with high precision as well as predicted unambiguously 
in the SM. 

• That a M can be measured so precisely, is kind of a miracle and possible only 
due to the specific properties of the muon. Due to the parity violating weak 
(V-A) interaction property, muons can easily be polarized and perfectly 
transport polarization information to the electrons produced in their decay. 

• There exists a magic energy ("magic 7") at which equations of motion 
take a particularly simple form. Miraculously, this energy is so high (3.1 
GeV) that the \i lives 30 times longer than in its rest frame! 



2 Even more promising would be a measurement of a r with additional enhancement 
(m T /m M ) 2 ~ 283. However, the much shorter lifetime of the r lepton (t t /t m ~ 
1.3 x 10~ 7 ) makes this measurement impossible at present. 
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□ /x's produced in pion decays are polarized 



\ / 

IP C TP 



□ polarized /i's decay producing electrons carrying the \i spin direction 




Fig. 1. Spin transfer properties in production and decay of the muons (P=parity, 
C=charge conjugation). fi~ [fi + ] is produced with positive [negative] helicity h = 
s -p/\p\, decay e~ [e + ] have negative [positive] helicity, respectively 

In fact only these highly energetic muons can by collected in a muon stor- 
age ring. At much lower energies muons could not be stored long enough to 
measure the precession precisely! 

Production and decay of the muons goes by the chain 

7T — > /J, + 

1 — ► e + v e + 

and the polarization "gymnastics" is illustrated in Fig.[TJ Note that the "max- 
imal" parity (P) violation means that the charged weak transition currents 
only couple to left-handed neutrinos u flL and right-handed antineutrinos v fJiR , 
in other words, parity violation is a direct consequence of the fact that the 
neutrinos v ilR and show no electromagnetic, weak and strong interaction 
in nature! as if they were non-existent. 



2 The BNL Muon g — 2 Experiment 



After the proposal of parity violation in weak transitions by Lee and Yang in 
1957, it immediately was realized that muons produced in weak decays of the 



6 The Muon g - 2 



pion (tt + — > neutrino) should be longitudinally polarized. In addition, 

the decay positron of the muon — > e + + 2 neutrinos) could indicate the 
muon spin direction. This was confirmed by Garwin, Lederman and Wein- 
rich 0] and Friedman and Telegdi [5jf|. The first of the two papers for the first 
time determined g^ = 2.00 within 10% by applying the muon spin precession 
principle. Now the road was free to seriously think about the experimental 
investigation of a M . 

The first measurement of {g^ — 2)/2 was performed at Columbia in 1960 [6] 
with a result a M = 0.00122(8) at a precsision of about 5%. Soon later in 1961, 
at the CERN cyclotron (1958-1962) the first precision determination became 
available [7]. Surprisingly, nothing special was observed within the 0.4% level 
of accuracy of the experiment. It was the first real evidence that the muon 
was just a heavy electron. In particular this meant that the muon is point-like 
and no extra short distance effects could be seen. This latter point of course 
is a matter of accuracy and the challenge to go further was evident. 

The idea of a muon storage ring was put forward next. A first one was 
successfully realized at CERN (1962-1968) [8]. It allowed to measure a M for 
both [i + and at the same machine. Results agreed well within errors and 
provided a precise verification of the CPT theorem for muons. An accuracy 
of 270 ppm was reached and an insignificant 1.7 a (1 a = 1 standard devi- 
ation) deviation from theory was found. Nevertheless the latter triggered a 
reconsideration of theory. It turned out that in the estimate of the three-loop 
0(a 3 ) QED contribution the leptonic "light-by- light scattering" part in the 
radiative corrections (dominated by the electron loop) was missing. Aldins et 
al. [9] then calculated this and after including it, perfect agreement between 
theory and experiment was obtained. 

The CERN muon g — 2 experiment was shut down end of 1976, while 
data analysis continued until 1979 |10) . Only a few years later, in 1984 the 
E821collaboration formed, with the aim to perform a new experiment at 
Brookhaven National Laboratory (BNL). Data taking was between 1998 and 
2001. The data analysis was completed in 2004. The E821 g~ 2 measurements 
achieved the remarkable precision of 0.54ppm [TTJ[T2], which is a 14-fold im- 
provement of the CERN result. The principle of the BNL muon g — 2 exper- 
iments involves the study of the orbital and spin motion of highly polarized 
muons in a magnetic storage ring. This method has been applied in the last 
CERN experiment already. The key improvements of the BNL experiment 
include the very high intensity of the primary proton beam from the Alter- 
nating Gradient Synchrotron (AGS), the injection of muons instead of pions 
into the storage ring, and a superferric storage ring magnet. The protons hit a 
target and produce pions. The pions are unstable and decay into muons plus 
a neutrino where the muons carry spin and thus a magnetic moment which 



3 The latter reference for the first time points out that P and C are violated simul- 
taneously, in fact P is maximally violated while CP is to very good approximation 
conserved in this decay (see Fig. [1} . 
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Protons Pions Polarized Muons 

from AGS H p=3.1 GeV/c Inflector -^Injection Point 




Fig. 2. The schematics of muon injection and storage in the g — 2 ring 

is directed along the direction of the flight axis. The longitudinally polarized 
muons from pion decay are then injected into a uniform magnetic field B 
where they travel in a circle (see Fig. [2]) . 




actual precession x 2 

Fig. 3. Spin precession in the g — 2 ring (~ 12° /circle) 

When polarized muons travel on a circular orbit in a constant magnetic 
field, as illustrated in Fig. [3l then a M is responsible for the Larmor preces- 
sion of the direction of the spin of the muon, characterized by the angular 
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10 20 30 40 50 60 70 80 



Time (ns) 

Fig. 4. Decay of fi + and detection of the emitted e + (PMT=Photomultiplier) 



fequency u) a . At the magic energy of about 
proportional to a^: 



3.1 GeV, the latter is directly 



,B 



1 



r 



- 1 



fix E 



E~3.1GcV 



at "magic 7" 



(9) 



Electric quadrupole fields E are needed for focusing the beam and they affect 
the precession frequency in general. 7 = E/m^ = 1/ yl — (3 2 is the relativistic 
Lorentz factor with (3 = v/c the velocity of the muon in units of the speed 
of light c. The magic energy E mag = 7 mag m /J is the energy E for which 
2 1 _ 1 = a a- The existence of a solution is due to the fact that a„ is a 
positive constant in competition with an energy dependent factor of opposite 
sign (as 7 > 1). The second miracle, which is crucial for the feasibility of the 
experiment, is the fact that 7mag = + a^)/a^ ~ 29.378 is large enough 
to provide the time dilatation factor for the unstable muon boosting the life 
time r M ~ 2.197 x 10 



717, ~ 6.454 x 10 sec, which allows 



to be stored in a ring of reasonable 



SeC tO Tin flight = 

the muons, traveling at v/c = 0.99942 • 
size (diameter ~ 14 m). 

This provided the basic setup for the g — 2 experiments at the muon storage 
rings at CERN and at BNL. The oscillation frequency u) a can be measured 
very precisely. Also the precise tuning to the magic energy is not the major 
problem. The most serious challenge is to manufacture a precisely known con- 
stant magnetic field B, as the latter directly enters the experimental extraction 
of ([9]) . Of course one also needs high enough statistics to get sharp values 
for the oscillation frequency. The basic principle of the measurement of a p is 
a measurement of the "anomalous" frequency difference aj a = \<^a\ — lu s — lu c , 
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20 40 60 80 100 

Time modulo 100ns [\is] 

Fig. 5. Distribution of counts versus time for the 3.6 billion decays in the 2001 
negative muon data-taking period. Courtesy of the E821 collaboration [11] 



where u> s = g^ (eh/2m fJ _) B/h = g M /2 • e/m M B is the muon spin-flip preces- 
sion frequency in the applied magnetic field and u> c = e/m^ B is the muon 
cyclotron frequency. The principle of measuring ui a is indicated in Fig. H and 
an example of a measured count spectrum is shown in Fig. [5l Instead of elim- 
inating the magnetic field by measuring u> c , B is determined from proton 
nuclear-magnetic-resonance (NMR) measurements. This procedure requires 
the value of /i^/ fi p to extract from the data. Fortunately, a high precision 
value for this ratio is available from the measurement of the hyperfine struc- 
ture in muonium. One obtains 

a » = I R \ 5 ' ( 10 ) 

where R — uj a /0j p and uj p = (e/m^c)(B) is the free-proton NMR frequency 
corresponding to the average magnetic field, seen by the muons in their or- 
bits in the storage ring. We mention that for the electron a Penning trap is 
employed to measure a e rather than a storage ring. The B field in this case 
can be eliminated via a measurement of the cyclotron frequency. The BNL 
g-2 muon storage ring is shown in Fig. [6l 

Since the spin precession frequency can be measured very well, the preci- 
sion at which g — 2 can be measured is essentially determined by the possi- 
bility to manufacture a constant homogeneous magnetic field B. Important 
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Fig. 6. The Brookhaven National Laboratory muon storage ring. The ring has a 
radius of 7.112 meters, the aperture of the beam pipe is 90 mm, the field is 1.45 
Tesla and the momentum of the muon is p M = 3.094 GeV/c. Picture taken from 
the Muon g — 2 Collaboration Web Page http://www.g-2.bnl.gov/ (Courtesy of 
Brookhaven National Laboratory) 



but easier to achieve is the tuning to the magic energy. The outcome of the 
experiment will be discussed later. 



3 QED Prediction of a e and the Determination of at 

The anomalous magnetic moment ag is a dimensionless quantity, just a num- 
ber, and corresponds to an effective tensor interaction term 

which in an external magnetic field at low energy takes the well known form 
of a magnetic energy (up to a sign) 



5£™ ^ _ Hm „ (jB (12) 

Such a term, if present in the fundamental Lagrangian, would spoil renormal- 
izability of the theory and contribute to Fyi{q 2 ) at the tree level. In addition, 
it is not SU(2)l gauge invariant, because gauge invariance only allows mini- 
mal couplings via a covariant derivative, i.e., vector and/or axial- vector terms. 
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The emergence of an anomalous magnetic moment term in the SM is a con- 
sequence of the symmetry breaking by the Higgs mechanism, which provides 
the mass to the physical particles and allows for helicity flip processes like 
the anomalous magnetic moment transitions. In any renormalizable theory 
the anomalous magnetic moment term must vanish at tree level. This means 
that there is no free adjustable parameter associated with it. It is a finite 
prediction of the theory. 

The reason why it is so interesting to have such a precise measurement of 
a e or a M , of course, is that it can be calculated with comparable accuracy in 
theory by a perturbative expansion in a of the form 

JV 

a t ~Y. A{2n) ( a l*) n > ( 13 ) 

n=l 

with up to TV = 5 terms under consideration at present. The experimental 
precision of a e (0.66 ppb) requires the knowledge of the coefficients with ac- 
curacies SA^ ~ 1 x 1CT 7 , SA^ - 6 x 1CT 5 , SA^ - 2 x 1CT 2 and 5A^ ~ 10. 
The expansion is an expansion in the number N of closed loops of the 
contributing Feynman diagrams. 

The recent new determination of a e |13j allows for a very precise determi- 
nation of the fine structure constant [HI Q~5] 

a _1 (a e ) = 137.035999069(96) [OJOppb] , (14) 

which we will use in the evaluation of . 

At two and more loops results depend on lepton mass ratios. For the eval- 
uation of these contributions precise values for the lepton masses are needed. 
We will use the following values for the muon-electron mass ratio, the muon 
and the tau mass [16j HZ] 

m^/me = 206.768 2838 (54) , m M /m r = 0.059 4592 (97) 
m e = 0.510 9989 918(44)MeV , m p = 105.658 3692 (94)MeV (15) 
m T = 1776.99 (29)MeV . 

The leading contributions to ai can be calculated in QED. With increasing 
precision higher and higher terms become relevant. At present, 4-loops are 
indispensable and strong interaction effects like hadronic vacuum polarization 
(vap) or hadronic light-by-light scattering (lbl) as well as weak effects have 
to be considered. Typically, analytic results for higher order terms may be 
expressed in terms of the Riemann zeta function 

CW =Ei (16) 

fc=i 

and of the poly-logarithmic integrals 
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Li n (ac) 



i-iy 



2)! 



\n n - 2 (t) ln(l-to) 
t 



dt = J2 f- 



fc=l 



(17) 



We first discuss the universal contributions in "one flavor QED", with 
one type of lepton lines only. At leading order one has 
• one 1-loop diagram 




a, = a,, 



2tt 



Schwinger 48 



giving the result mentioned before. 

• At 2-loops 7 diagrams with only (-type fermion lines 










which contribute a term 

(4) 

a) 



197 
144 



7T 

12 



In 2 



rC(3) 



© 2 . 



(18) 



obtained independently by Peterman [18j and Sommerfield [19j in 1957. 

• At 3-loops, with one type of fermion lines only, the 72 diagrams of Fig. [7] 
contribute. Most remarkably, after about 25 years of hard work, Laporta and 
Remiddi in 1996 [3U] managed to give a complete analytic result (see also [2Tj ) 



,( 6 ) 



28259 17101 



298 



5184 
1 , 



24 



810 



In 2 2 )■ - 



ir 2 In 2 



139 
l8~ 



C(3) 



100 
IT 



Li 4 (I) + ^ln 4 2 
y 2' 24 



239 
2160 7 



(f)' 



(19) 



This result was confirming Kinoshita's earlier numerical evaluation [22 . 

The big advantage of the analytic result is that it allows a numerical 
evaluation at any desired precision. The direct numerical evaluation of the 
multidimensional Feynman integrals by Monte Carlo methods is always of 
limited precision and an improvement is always very expensive in computing 
power. 

• At 4-loops 891 diagrams contribute to the universal term. Their evalua- 
tion is possible by numerical integration and has been performed in a heroic 
effort by Kinoshita [53] (reviewed in [H]), and was updated recently by Ki- 
noshita and his collaborators (2002/2005/2007) [2S1[T5]. 

The largest uncertainty comes from 518 diagrams without fermion loops 

(8) 

contributing to the universal term A\ . Completely unknown is the universal 
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Fig. 7. The universal third order contribution to a M . All fermion loops here are 
muon-loops (first 22 diagrams). All non-universal contributions follow by replacing 
at least one muon in a closed loop by some other fermion 



five-loop term A^ 10 ', which is leading for a e . An estimation discussed in |39j 
suggests that the 5-loop coefficient has at most a magnitude of 3.8. We adopt 
this estimate and take into account A^ 10 -* = 0.0(3.8) (as in [25]). 
Collecting the universal terms we have 

a} ni = 0.5 (-) - 0.32847896557919378 ... (-) 2 

+ 1.181241456587... (-J - 1.9144(35) [-) +0.0(3.8) [-) 
= 0.001 159652 176 42(81) (10) (26) [86] • ■ • (20) 

for the one-flavor QED contribution. The three errors are: the error of a, given 
in (fT4"]) , the numerical uncertainty of the a 4 coefficient and the estimated size 
of the missing higher order terms, respectively. 

At two loops and higher, internal fermion-loops show up, where the flavor 
of the internal fermion differs form the one of the external lepton, in general. 
As all fermions have different masses, the fermion-loops give rise to mass 
dependent effects, which were calculated at two-loops in [26l [27] (see also [28l 
OH [301 [31]), and at three-loops in [3l [33] [M E3 [33 [37] . 

The leading mass dependent effects come from photon vacuum polariza- 
tion, which leads to charge screening. Including a factor e and considering 
the renormalized photon propagator (wave function renormalization factor 
Z^) we have 
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i e 2 D^g) = - 2 (f" e n/ fl^ + gauge terms (21) 

which in effect means that the charge has to be replaced by an energy- 
momentum scale dependent running charge 

^^ = T^F Y (22) 

The wave function renormalization factor Z 7 is fixed by the condition that 
as q 2 — » one obtains the classical charge (charge renormalization in the 
Thomson limit). Thus the renormalized charge is 

eW <^= i + w( ;U;(o)) (23) 

where in perturbation theory the lowest order diagram which contributes to 

n'^q 2 ) is 

7 1 



L 

and describes the virtual creation and re-absorption of fermion pairs 7* — > 
e + e~ , , t + t~ , mm, cM, • • • (had) — > 7* . 

2 

In terms of the fine structure constant a = |- Eq. (|23j) reads 

"(g 2 )= l_^ a(g2) ! ^(g 2 ) = -Re (7i;(g 2 )- 77^(0)) . (24) 

The various contributions to the shift in the fine structure constant come 
from the leptons (lep = e, [i and r), the 5 light quarks (u, b, s, c, and b) 
and/or the corresponding hadrons (had). The top quark is too heavy to give 
a relevant contribution. The hadronic contributions will be considered later. 
The running of a is governed by the renormalization group (RG). In the 
context of g — 2 calculations, the use of RG methods has been advocated 
in [30]. In fact, the enhanced short-distance logarithms may be obtained by 
the substitution a — ► a(m^) = a (1 + 3^ In ^f- + • • • ) in a lower order result 
(see the following example). 

Typical contributions are the following: 
— LIGHT internal masses give rise to log's of mass ratios which become sin- 
gular in the light mass to zero limit (logarithmically enhanced corrections) 



/6\ 



3 m e 36 y J 



— HEAVY internal masses decouple, i.e., they give no effect in the heavy 
mass to infinity limit 
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A = [M^) 2+o 0l ln ^) 

New physics contributions from states which are too heavy to be produced 
at present accelerator energies typically give this kind of contribution. Even 
so et M is 786 times less precise than a e it is still 54 times more sensitive to new 
physics (NP). 

Corrections due to internal e, /i- and r-loops are different for a e , and a T . 
For reasons of comparison and because of its role in the precise determination 
of a we briefly consider a e first. The result is of the form 

a QED = a uni + fle(/i) + fle(r) + t) (25) 

witrfl [13 [351 \M 

-J - 7.373 941 64(29) x 1CT 6 (-J 




a e (r) = 1.83763(60) x 10 9 (^— J - 6.5819(19) x 10 8 (^— J 



a e (/x,r) = 0.190945(62) x 10~ 12 ' 
The QED part thus may be summarized in the prediction 
oed = _ 0.328 478 444 002 90(60) f-V 

27T \7T/ 

+ 1.181234016828(19) ^- J - 1.9144(35) {- J +0.0(3.8) (-J . (26) 

The hadronic and weak contributions to a e are small : cig ad — 1.67(3) x 10 -12 
and a™ ak = 0.036 x 10 -12 , respectively. The hadronic contribution now just 
starts to be significant, however, unlike in ajj ad for the muon, a^ ad is known 
with sufficient accuracy and is not the limiting factor here. The theory error 
is dominated by the missing 5-loop QED term. As a result a e essentially only 



The order a 3 terms are given by two parts which cancel partly 

4 S1 (m«/%) = -2.17684015(11) x 10~ 5 | M _ vap + 1.439445989(77) x 10" 5 | M _ lbl 
A ( 2 6 \m e /m T ) = -1.16723(36) x 10~ 7 | T _ vap + 0.50905(17) x 10" 7 | t _ 1w . 

The errors are due to the uncertainties in the mass ratios. They are negligible 
in comparison with the other errors, "vap" denotes vacuum polarization type 
contributions [35] and "lbl" light-by-light scattering type ones [36] (the first 6 
diagrams of Fig. [7] with an e- or r- loop). 
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depends on perturbative QED, while hadronic, weak and new physics (NP) 
contributions are suppressed by (m e /M) 2 , where M is a weak, hadronic or 
new physics scale. As a consequence a e at this level of accuracy is theoretically 
well under control (almost a pure QED object) and therefore is an excellent 
observable for extracting oqed based on the SM prediction 

af M = a^ ED [Eq. ®] + 1.706(30) x 10~ 12 (hadronic & weak) . (27) 

We now compare this result with the very recent extraordinary precise 
measurement of the electron anomalous magnetic momen10 [13j 

a^ xp = 0.001 159 652180 85(76) (28) 

which yields 

a _1 (a e ) = 137.035999069(90)(12)(30)(3) , 

which is the value (fT4|) [14j[15] we use in calculating a M . The first error is the 
experimental one of a™ p , the second and third are the numerical uncertainties 
of the a 4 and a 5 terms, respectively. The last one is the hadronic uncertainty, 
which is completely negligible. Note that the largest theoretical uncertainty 
comes from the almost completely missing information concerning the 5-loop 
contribution. This is now the by far most precise determination of a and we 
will use it throughout in the calculation of a M , below. 
The best determinations of a which do not depend on a e are [401 141] 

a _1 (Cs06) = 137.03600000(110)[8.0ppb] , 
a _1 (Rb06) = 137.03599884(091) [6.7 ppb] , 

less precise by about a factor ten. «(Cs06) is determined from a measurement of 
h/Mcs via Cesium recoil measurements 40 , while a(Rb06) derives from the ratio 
h/Alub measured via Bloch oscillations of Rubidium atoms in an optical lattice [41] . 
These values should be used in theoretical predictions of a e . Using a(Cs06) we get 
a e = 0.00115965217298(930) and of* - a* he = 7.87(9.33) x 10~ 12 , with «(Rb06) 
the prediction reads a e = 0.00115965218279(769) and of - af* = -1.94(7.73) x 
10~ 12 in best agreement. The error of the prediction is completely dominated by 
the uncertainty coming from a(Cs06) and a(Rb06) such that an improvement of a 
by a factor 10 would allow a much more stringent test of QED (see also [141 115] ). If 
one assumes that |/\a^ ew p "y slcs ~ ml/ A 2 where A approximates the scale of "New 
Physics", then the agreement between a~ 1 (a e ) and « -1 (Rb06) currently probes 
A £j 0(250 GeV). To access the much more interesting A ~ 0(1 TeV) region also a 
bigger effort on the theory side would by necessary about the 0(a 4 ) and the 0(a 5 ) 
terms. 



5 The famous g e measurement from University of Washington (Dehmelt et al. 
1987) [38] found af p = 0.001 159 652 188 30(420) and recently has been improved 
by about a factor 6 in an experiment at Harvard University (Gabrielse et al. 
2006). The new central value shifted downward by 1.7 standard deviations. 
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4 Standard Model Prediction for a. 



4.1 QED Contribution 

The SM prediction of looks formally very similar to the one for a e , however, 
besides the common universal part, the mass dependent, the hadronic and 
the weak effects enter with very different weight and significance. The mass- 
dependent QED corrections follow from the universal set of diagrams (see 
e.g. Fig. [7] for the 3 loop case) by replacing the closed internal /x-loops by 
e- and/or r-loops. Typical contributions come from vacuum polarization or 
light-by-light scattering loops, like 




4 6) (lbl,e) 



-7T In h — TT 

3 m e 270 



3C(3) 



3 3 \m u m e 



The result is given by 

a M = a" m + a Al (m M /m e ) + a Al (m Al /m T ) + a M (m M /m e , m^/m T ) (29) 
witlfl [271 [351 EH [37] 

aJmJm e ) = 1.094 258 311 1 (84) (-) +22.868 380 02 (20) (- 

+ 132.682 3 (72) 

a M (TO M /m T ) = 7.8064 (25) x 10~ 5 (-) 2 + 36.051 (21) x 10~ 5 (-) 

/a\ 4 

+ 0.005 (3) (-) 



ttcK/me.mp/frv) = 52.766 (17) x 10"° ( - ] 0.037 5!>4 



except for the last term, which has been worked out as a series expansion in the 
mass ratios [42 , 43 , all contributions are known analytically in exact form [35l 



6 Again the order a 3 terms are given by two parts (see {T3}) 

Af(ra P /"" e ) = 20.94792489(16)| e _ lbl + 1.920 455130(33)1 



A< 6) (m M /m r ) = 0.002 142 83(69)| T _ 



0.001782 33(48) | T _ 



The errors are due to the uncertainties in the mass ratios. Note that the electron 
light- by-light scattering loop gives an unexpectedly large contribution [5]. 
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up to 3-loops. At 4-loops only a few terms are known analytically 45J. 
Again the relevant 4-loop contributions have been evaluated by numerical 
integration methods by Kinoshita and Nio [IB]. The 5-loop term has been 
estimated to be A^ 10) (m M /m e ) = 663(2) in [4711151115) . 

Our knowledge of the QED result for may be summarized by 

oed = at_ + Q 7g5 41Q(26) / «\ 2 
F 2n \7r/ 

+24.050 509 65(46) (^) 3 + 130.8105(85) + 663(20) (^) 5 . (30) 

Growing coefficients in the a/ir expansion reflect the presence of large In — ^ ~ 
5.3 terms coming from electron loops. In spite of the strongly growing expan- 
sion coefficients the convergence of the perturbation series is excellent 



# n of loops 


Q [(a/ir) n ] 


x 10 11 


1 


+0.5 


116140973.30 (0.08) 


2 


+0.765 857410(26) 


413217.62 (0.01) 


3 


+24.050 509 65(46) 


30141.90 (0.00) 


4 


+130.8105(85) 


380.81 (0.03) 


5 


+663.0(20.0) 


4.48 (0.14) 


tot 




116584718.11 (0.16) 



because a /it is a truly small expansion parameter. 

The different higher order QED contributions are collected in Tab. [TJ We 
thus arrive at a QED prediction of given by 



Table 1. QED contributions to a M in units 10 



term 


universal 


e-loops 


r-loops 


e&T-loops 


a« 
a< 6 > 
a W 
a' 10 ' 


-1.772 305 06 (0) 
0.014 804 20 (0) 

-0.000 055 73(10) 
0.000 000 00(26) 


5.904 060 07 (5) 
0.286 603 69 (0) 
0.003 862 56 (21) 
0.000 044 83(135) 


0.000 42120(13) 
0.000 004 52 (1) 
0.000 00015 (9) 

? 


0.000 006 61(0) 
0.000 00109(0) 

? 



g qed = H6 584 718.113(.082)(.014)(.025)(.137)[.162] x 10" 11 (31) 

where the first error is the uncertainty of a in (|14|) , the second one combines 
in quadrature the uncertainties due to the errors in the mass ratios, the third 
is due to the numerical uncertainty and the last stands for the missing 0(a 5 ) 
terms. With the new value of a[a e ] the combined error is dominated by our 
limited knowledge of the 5-loop term. 

7 Explicitly, the papers only present expansions in the mass ratios; some result have 
been extended in [37] and cross checked against the full analytic result in |44] , 



4.2 Weak Contributions 
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The electroweak SM is a non-Abelian gauge theory with gauge group 51/ (2) 
C/(l)y — > C/(1)qed, which is broken down to the electromagnetic Abelian 
subgroup C/(1)qed by the Higgs mechanism, which requires a scalar Higgs 
field H which receives a vacuum expectation value v. The latter fixes the 
experimentally well known Fermi constant G M = l/(v / 2w 2 ) and induces the 
masses of the heavy gauge bosons Mw and Mz as well as all fermion masses 
m / . Other physical constants which we will need later for evaluating the weak 
contributions are the Fermi constant 

G M = 1.16637(1) x 1(T 5 GeV~ 2 , (32) 

the weak mixing parameter 

sin 2 O w = 0.22276(56) (33) 

and the masses of the intermediate gauge bosons Z and W 

M z = 91.1876 ± 0.0021 GeV , M w = 80.392 ± 0.029 GeV . (34) 

For the not yet discovered SM Higgs particle the mass is constrained by LEP 
data to the range 

114 GeV <m H < 200 GeV (at 96% CL) . (35) 

The weak interaction contributions to a M are due to the exchange of the 
heavy gauge bosons, the charged and the neutral Z, which mixes with 
the photon via a rotation by the weak mixing angle Ow and which defines 
the weak mixing parameter sin 2 Ow = 1 — M w /M z . What is most interesting 
is the occurrence of the first diagram of Fig. [51 which exhibits a non-Abelian 
triple gauge vertex and the corresponding contribution provides a test of the 
Yang-Mills structure involved. It is of course not surprising that the photon 
couples to the charged W boson the way it is dictated by electromagnetic 
gauge invariance. The gauge boson contributions up to negligible terms of 




Fig. 8. The leading weak contributions to at, diagrams in the physical unitary gauge 

m 2 

order 0( Mi " ) are given by [SU] 
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,(2) EW 



(W) = 



VZG^ml 10 



167T 



, — ~ +388.70(0) x 10" 
z 3 



EW, 



a)7'"'\Z) 



y/2Gpm* (-1 + 4 sin 2 O w f - 5 



-193.88(2) x 10 



-li 



16tt 2 3 
while the diagram with the Higgs exchange, for tojj 3> to„, yields 

2) EW / 



'(H) 



V2G M m 2 , to 2 



4tt 2 



— £- In — £- + ■ • • < 5 x 10~ 14 for m H > 114 GeV 



""H H 

Employing the SM parameters given in ((32)) and (|33|) we obtain 



a 



(2) EW 



(194.82 ±0.02) x 10 



-li 



(36) 



The error comes from the uncertainty in sin 2 Ow given above. 

The electroweak two-loop corrections have to be taken into account as 
well. In fact triangle fermion-loops may give rise to unexpectedly large radia- 
tive corrections. The diagrams which yield the leading corrections are those 
including a VVA triangular fermion-loop (VVA ^ while VVV — 0) associ- 
ated with a Z boson exchange 

7 




which exhibits a parity violating axial coupling (A). A fcrmion of flavor / 
yields a contribution 



(4) EW 



([/]) 



m 2 , a 



16tH 



2T 3f N cf Q 



M 2 

31n— f + C f 

m 2 , 



(37) 



where T3 / is the 3rd component of the weak isospin, Q / the charge and N c f the 
color factor, 1 for leptons, 3 for quarks. The mass to// is if to/ < and 
m/ if TO/ > to m , and C e = 5/2, C'^ = 11/6 — 8/9 7r 2 , C T = —6 [51] . However, 
in the SM the consideration of individual fermions makes no sense and a 
separation of quarks and leptons is not possible. Mathematical consistency 
of the SM requires complete VVA anomaly cancellation between leptons and 
quarks, and actually J^f ^cfQ 2 jT^f — holds for each of the 3 known lepton- 
quark families separately. Treating, in a first step, the quarks like free fermions 
(quark parton model QPM) the first two families yield (using m u — m c i — 
300 MeV , m s = 500 MeV , m c = 1.5GeV) 



(4) EW 



e, u, d 
H,c,s 



)qpm - 



ml a 



16tH 



7T 



In 



19 2 2 

m ti m d m s 



49 

y 



87T 2 

9 



16tt 



tti 2 , a 

5-& - X 32.0(?) 



-8.65(?) x 10" 



(38) 
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which demonstrates that the leading large logs ~ In Mz have canceled |52) . 
as it should be. However, the quark masses which appear here are ill-defined 
constituent quark masses, which can hardly account reliably for the strong 
interaction effects, therefore the question marks in place of the errors. 

In fact, low energy QCD is characterized in the chiral limit of massless 
light quarks u, d, s, by spontaneous chiral symmetry breaking (SxSB) of the 
chiral group SU(3)v ® SU(3)a, which in particular implies the existence of 
the pseudoscalar octet of pions and kaons as Goldstone bosons. The light 
quark condensates are essential features in this situation and lead to non- 
perturbative effects completely absent in a perturbative approach. Thus low 
energy QCD effects are intrinsically non-perturbative and controlled by chi- 
ral perturbation theory (CHPT), the systematic QCD low energy expansion, 
which accounts for the S^SB and the chiral symmetry breaking by quark 
masses in a systematic manner. The low energy effective theory describing 
the hadronic contributions related to the light quarks u, d, s requires the cal- 
culation of the diagrams of the type shown in Fig. [9] The leading effect for 




(a) [L.D.] (b) [L.D.] (c) [S.D.] 

Fig. 9. The two leading CHPT diagrams (L.D.) and the QPM diagram (S.D.). The 
charged pion loop is sub-leading and is discarded. Diagrams with permuted 7 <-> Z 
on the /i-line have to be included 



the 1st plus 2nd family takes the form 



(4) EW 



e, u, d 
/!, c, s 



V2G, 



mf, a 



;chpt 



16tt 2 



_14 1q M| 
3 mi 



41„M 
1 - 



i)ir 



35 

y 



— TT 

9 



16tt 2 



x 26.2(5) ~ -7.09(13) x 10 



-11 



(39) 



The error comes from varying the cut-off Ma between l GeV and 2 GeV. 
Below l GeV CHPT can be trusted above 2 GeV we can trust pQCD. For- 
tunately the result is not very sensitive to the choice of the cut-off. For more 
sophisticated analyses we refer to [52j [53l [54] which was corrected and re- 
fined in [55, 56]. Thereby, a new kind of non-renormalization theorems played 
a key role (57] EH [59]. Including subleading effects yields —6.7 x 10~ n for 
the first two families. The 3rd family of fermions including the heavy top 
quark can be treated in perturbation theory and has been worked out to 
be —8.2 x 10~ n in [60]. Subleading fermion loops contribute —5.3 x 10~ n . 
There are many more diagrams contributing, in particular the calculation 
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of the bosonic contributions (1678 diagrams) is a formidable task and has 
been performed 1996 by Czarnecki, Krause and Marciano as an expansion in 
(m^/My) 2 and (My /tuh) 2 [SI]- Later complete calculations, valid also for 
lighter Higgs masses, were performed |62( 163] . which confirmed the previous 
result —22.3 x 10 -11 . The 3rd family of fermions including the heavy top quark 
can be treated in perturbation theory and has been worked out in [60] . 
The complete weak contribution may be summarized by |56j 



with errors from triangle quark loops and from variation of the Higgs mass 
in the range win = 150l^° GeV. The 3-loop effect has been estimated to be 
negligible [55ll5rJ| . 

4.3 Hadronic Contributions 

So far when we were talking about fermion loops we only considered the lep- 
ton loops. Besides the leptons also the strongly interacting quarks have to be 
taken into accounl[f|. The problem is that strong interactions at low energy 
are non-perturbative and straight forward first principle calculations become 

very difficult and often impossible. 

Fortunately the leading hadronic effects are vacuum polarization type cor- 
rections (see ([2"3")) ). which can be safely evaluated by exploiting causality (an- 
alyticity) and unitarity (optical theorem) together with experimental low en- 
ergy data. In fact vacuum polarization effects may be calculated using the 
master formula 



which replaces a free photon propagator by a dressed one, and where the 
imaginary part of the photon self-energy function iT 7 (s) is determined via the 
optical theorem by the total cross-section of hadron production in electron- 
positron annihilation: 

8 The theory of strong interactions is Quantum Chromodynamics (QCD) |64] , The 
strongly interacting particles, the hadrons, are made out of quarks and/or anti- 
quarks, which interact via an octet of gluons according to the non-Abelian SU(3) C 
gauge theory. The gauged internal degrees of freedom are named color. Quarks 
are flavored and labeled as up (it), down (d), strange (s), charm (c), bottom (b) 
and top (t). Each of the flavored quarks exists in N c = 3 colors (red, green, 
blue). All hadrons are color neutral bound states (confinement). This means that 
QCD is intrinsically non-perturbative. However, QCD also has the property of 
asymptotic freedom [BS], which implies that perturbation theory starts to work 
at higher energies, where the quark structure appears resolved as in deep inelastic 
electron-proton scattering, for example. 




(154±l[had]±2[m H ,m t ,3-loop]) x 10 



li 



(40) 




(41) 



Essentials of the Muon g — 2 23 
47r 2 a 1 

cr(s) e+e -^ 7 .^ hadrons = Im 77 7 (s) . (42) 

S 7T 

The leading hadronic contribution is represented by the diagram Fig. 1 1 01 




Fig. 10. The leading order (LO) hadronic vacuum polarization diagram 



which corresponds to a contribution a ™ asslvc ~f = K(s) of the lowest order 
diagram with the photon replaced by a "massive photon" of mass ^fs, and 
convoluted according to (|4T|) . It yields the dispersion integral 



a f°° ds 1 , . . . , , f 1 , a; 2 (l- x) , s 

<V = ~ / ImlT 7 (s)if(s) , iT(s) = / dx^ v y (43) 



7T 



II 



As a result the leading non-perturbative hadronic contributions a^ ad can be 

obtained in terms of R-y(s) = cr'°'(e + e~ — > 7* — > hadrons) / 4 ^° data via the 
dispersion integral: 



had 



(^) 2 ( J\ s <^ s Kj ds R ^ c y^ (44) 



The rescaled kernel function K(s) = is/m 2 ^ K{s) is a smooth bounded func- 
tion, increasing from 0.63... at s = Am 2 , to 1 as s — ► 00. The 1/s 2 enhancement 
at low energy implies that the p — > 7r + 7r~ resonance is dominating the dis- 
persion integral (~ 75 %). Data can be used up to energies where 7 — Z 
mixing comes into play at about 40 GeV. However, by the virtue of asymp- 
totic freedom, perturbative Quantum Chromodynamics (pQCD) becomes the 
more reliable the higher the energy and in fact may be used safely in re- 
gions away from the flavor thresholds where the non-perturbative resonances 
show up: p, u>, </>, the J /if) series and the T series. We thus use perturbative 
QCD [551 [57] from 5.2 to 9.46 GeV and for the high energy tail above 13 GeV, 
as recommended in [55J [571 EE] • 

Hadronic cross section measurements e + e~ — > hadrons at electron-positron 
storage rings started in the early 1960's and continued up to date. Since our 
analysis [69] in 1995 data from MD1 [70], BES-II [71] and from CMD-2 [72] 
have lead to a substantial reduction in the hadronic uncertainties on a^ ad . 
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More recently, KLOE [73], SND [71] and CMD-2 [75] published new measure- 
ments in the region below 1.4 GeV. My up-to-date evaluation of the leading 
order hadronic VP yields [75] 

a had(i) = ( 6g2 A ± 5>6 j x 10 -io _ (45) 
Some other recent evaluations are collected in Tab. [2] Differences in errors 



Table 2. Some recent evaluations of a 



had(l) 



a^ ad(1) X 10 10 


data 




Ref. 


696.3(7.2] 


e + e~ 




[77] 


711.0(5.8] 


e + e~ 


+ T 


[77] 


694.8(8.6] 


e + e~ 




EH] 


684.6(6.4] 


e + e~ 


TH 


79 


699.6(8.9] 


e + e~ 




80 


692.4(6.4] 


e + e~ 




ED 


693.5(5.9] 


e + e~ 




82 


701.8(5.8] 


e + e~ 


+ T 


82 


690.9(4.4] 


e + e~ 




[Ha 


689.4(4.6] 


e + e~ 




84 


692.1(5.6] 


e + e~ 







Davier et al. (03) (e + e ) 
Davier et al. (03) (e + e~,r) 
Ghozzi, Jegerlehner (03) (e + e ) 
Narison (03) (e + e") 
Ezhela et al. (03) (e + e") 
Hagiwara et al. (03) (e + e~)(ird) 
Troconiz, Yndurain (04) (e + e~) 
Troconiz, Yndurain (04) (e + e~,r) 
Davier et al. (06) (e + e") 
Hagiwara et al. (06) (e + e~)(ircl) 
Jegerlehner (06) (e + e~) 



680.0 700.0 720.0 



come about mainly by utilizing more "theory-driven" concepts : use of selected 
data sets only, extended use of perturbative QCD in place of data [assuming 
local duality] , sum rule methods and low energy effective methods (85 . Only 
the last three (**) results include the most recent data from SND, CMD-2, 
and BaBai0. 

In principle, the 1=1 iso-vector part of e + e~ — * hadrons can be ob- 
tained in an alternative way by using the precise vector spectral functions 
from hadronic r-decays r — > v T + hadrons which are related by an isospin 
rotation [86j . After isospin violating corrections, due to photon radiation and 
the mass splitting rrid — rn u ^ 0, have been applied, there remains an unex- 
pectedly large discrepancy between the e + e~- and the r-based determinations 
of a M [77] , as may be seen in Table [2] Possible explanations are so far unac- 
counted isospin breaking [78] or experimental problems with the data. Since 

9 The analysis [84] does not include exclusive data in a range from 1.43 to 2 GeV; 
therefore also the new BaBar data are not included in that range. It also should 
be noted that CMD-2 and SND are not fully independent measurements; data 
are taken at the same machine and with the same radiative correction program. 
The radiative corrections play a crucial role at the present level of accuracy, and 
common errors have to be added linearly. In I77!,!83_ pQCD is used in the extended 
ranges 1.8 - 3.7 GeV and above 5.0 GeV; furthermore |83j excludes the KLOE 
data. 
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the e + e~-data are more directly related to what is required in the dispersion 
integral, one usually advocates to use the e + e~ data only. 

At order 0(a 3 ) diagrams of the type shown in Fig. [TT] have to be cal- 
culated, where the first diagram stands for a class of higher order hadronic 
contributions obtained if one replaces in any of the first 6 two-loop diagrams 
on p. [12] one internal photon line by a dressed one. The relevant kernels for 




Fig. 11. Higher order (HO) vacuum polarization contributions 

the corresponding dispersion integrals have been calculated analytically in [87\ 
and appropriate series expansions were given in [88] (for earlier estimates see 
[HS1 ED] ). Based on my recent compilation of the e + e~ data [76] I obtain 

a had( 2 ) = (_ioo.3±2.2) x HT 11 , (46) 

in accord with previous/other evaluations [90 | l88 [ l86" [ l81~ l [84] . 

Much more serious problems with non-perturbative hadronic effect we en- 
counter with the hadronic light-by-light (LbL) contribution at 0(a 3 ) depicted 
in Fig. 1121 Experimentally, we know that 77 — > hadrons — > 77 is dominated 




Fig. 12. Hadronic light- by- light scattering in g — 2 



by the hadrons ir°, rj, 77', • • • , i.e., single pseudoscalar meson spikes [91], and 
that 7T° — * 77 etc. is governed by the parity odd Wess-Zumino-Witten (WZW) 
effective Lagrangian 

£(4) = __^c F ^ A P d * n + . . . (47) 

12 vr/o 

which reproduces the Adler-Bell-Jackiw triangle anomaly and which helps in 
estimating the leading hadronic LbL contribution, fo denotes the pion decay 
constant /„• in the chiral limit of massless light quarks. Again, in a low energy 
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(a) [L.D.] (b) [L.D.] (c) [S.D.] 

Fig. 13. Leading hadronic light-by-light scattering diagrams: the two leading 
CHPT diagrams (L.D.) and the QPM diagram (S.D.). The charged pion loop is 
sub-leading only, actually. Diagrams with permuted 7's on the p-line have to be 
included. 7-hadron/quark vertices at q 2 7^ are dressed (VMD) 



effective description, the quasi Goldstone bosons, the pions and kaons play an 
important role, and the relevant diagrams are displayed in Fig[T3l 

However, as we know from the hadronic VP discussion, the p meson is 
expected to play an important role in the game. It looks natural to apply a 
vector-meson dominance (VMD) like model. Electromagnetic interactions of 
pions treated as point-particles would be descried by scalar QED in a first step. 
However, due to hadronic interactions the photon mixes with hadronic vector- 
mesons like the p°. The naive VMD model attempts to take into account this 
hadronic dressing by replacing the photon propagator as 

iff"" m l , 
q 2 q 2 — TOp 

where the ellipses stand for the gauge terms. The main effect is that it provides 
a damping at high energies with the p mass as an effective cut-off (physical 
version of a Pauli-Villars cut-off). However, the naive VMD model is not com- 
patible with chiral symmetry. The way out is the Resonance Lagrangian Ap- 
proach (RLA) 92J , an extended version of CHPT which incorporates vector- 
mesons in accordance with the basic symmetries. The Hidden Local Symmetry 
(HLS) model and the Extended Nambu-Jona-Lasinio (ENJL) model 
are alternative versions of RLA, which are basically equivalent [95 , for what 
concerns this application. 

Based on such effective field theory (EFT) models, two major efforts in 
evaluating the full a^ bL contribution were made by Hayakawa, Kinoshita 
and Sanda (HKS 1995) [H], Bijnens, Pallante and Prades (BPP 1995) [57] 
and Hayakawa and Kinoshita (HK 1998) [5H] (see also Kinoshita, Nizic and 
Okamoto (KNO 1985) [90]). Although the details of the calculations are quite 
different, which results in a different splitting of various contributions, the re- 
sults are in good agreement and essentially given by the 7r°-pole contribution, 
which was taken with the wrong sign, however. In order to eliminate the cut- 
off dependence in separating L.D. and S.D. physics, more recently it became 
favorable to use quark-hadron duality, as it holds in the large N c limit of 
QCD, for modeling of the hadronic amplitudes [9H]. The infinite series of nar- 



1.9 



1.9 



q 2 — m 2 
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row vector states known to show up in the large N c limit is then approximated 
by a suitable lowest meson dominance (LMD+V) ansatz |100j , assumed to be 
saturated by known low lying physical states of appropriate quantum num- 
bers. This approach was adopted in a reanalysis by Knecht and Nyffeler (KN 
2001) [TOll EH HH3 HD3 in 2001, in which they discovered a sign mistake 
in the dominant 7r°, 77, 77' exchange contribution, which changed the central 
value by +167 x 10~ n , a 2.8 a shift, and which reduces a larger discrepancy 
between theory and experiment. More recently Melnikov and Vainshtein (MV 
2004) [105] found additional problems in previous calculations, this time in 
the short distance constraints (QCD/OPE) used in matching the high energy 
behavior of the effective models used for the ir°, 77, rf exchange contribution. 

ThcP^I most important pion-pole term is of the form (p is the muon mo- 
mentum, qi (i — 1,2,3) are the virtual photon momenta, two of which are 
chosen as loop integration variables) [101] 



< -° = -e 6 



dV_ d^2_ 1 

(2^) 4 (2^) 4 q\ql(q x + q 2 ) 2 [(p + qi ) 2 - m 2 ][(p - q 2 ) 2 - m< 

I „1 „2 „2\ t I „2 „2 n\ 



T 7r*7*7* (?2' 9l ' ?!) Tr*l'lil 2 ■ 92' 0) „ / n 

— 2 — 2 1 ■ T 1 (q 1 ,q 2 ;p) 

•? 7 7r*7*7* q%, q 2 ) ^•7*7(93, ?3, 0) . , 

2 2——. -L2{qi,q2,p) 



q$-m 2 .+ IE 



(48) 



where Xi (gi , #2 ; and ^((Zi) <?2jp) are known scalar kinematics factors and 
gf ) is the non-perturbative 7r°77 form factor (FF) whose off- 
shell form is essentially unknown in the integration range of (|48|l . 

A new quality of the problem encountered here is the fact that the inte- 
grand depends on 3 invariants q 2 , q 2 , q 2 q^ — —(qi + q 2 ). While hadronic VP 
correlators or the VVA triangle with an external zero momentum vertex only 
depends on a single invariant q 2 . In the latter case the invariant amplitudes 
(form factors) may be separated into a low energy part q 2 < A 2 (soft) where 
the low energy effective description applies and a high energy part q 2 > A 2 
(hard) where pQCD works. In multi-scale problems, however, there are mixed 
soft-hard regions where no answer is available in general, unless we have data 
to constrain the amplitudes in such regions. In our case, only the soft region 
9i ) 92 j 93 — ^ 2 an d the hard region q\,q\,q\ > A 2 are under control of ei- 
ther the low energy EFT and of pQCD, respectively. In the mixed soft-hard 
domains operator product expansions and/or soft versus hard factorization 
"theorems" a la Brodsky-Farrar [106j may help. Actually, one more approxi- 
mation is usually made: the pion-pole approximation ,i.e., the pion- momentum 
square (first argument of F) is set equal to m 2 , as the main contribution is ex- 
pected to come from the pole. Knecht and Nyffeler modeled T^i* (w 2 , q 2 , q 2 ) 
in the spirit of the large N c expansion as a "LMD+V" form factor: 



This paragraph cannot be more than a rough sketch of an ongoing discussion. 
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T , 2 2 2 s U llltill + 92 ) + Mgl + 92) 2 + ^gj^ + Mg| + gf) + ft7 

^,,,(m„Ji,fc)- 3 ( g f-Mf)(gf-M|)( g l-M?)(gl-M|) ' ^> 

with /i 7 = -(N c MfMi/Air 2 f 2 ), / w ~ 92.4 MeV. An important constraint 
comes from the pion-pole form factor ^ 7 » 7 (m 2 , — Q 2 , 0), which has been 
measured by CELLO [107] and CLEO |108j . Experiments are in fair agreement 
with the Brodsky-Lepage |109j form 

^ - Q2 ' 0) ^ - ^uirm^m (50) 

which interpolates between a l/Q 2 asymptotic behavior and the constraint 
from 7T° decay at Q 2 — 0. This behavior requires hi — 0. Identifying the 
resonances with Mi = M p = 769 MeV, M 2 = M p > = 1465 MeV, the phe- 
nomenological constraint fixes /15 = 6.93 GeV . h 2 will be fixed by later. 
As the previous analyses, Knecht and Nyffeler apply the above VMD type 
form factor on both ends of the pion line. In fact at the vertex attached to 
the external zero momentum photon, this type of pion-pole form factor can- 
not apply for kinematical reasons: when q^ xt = not J r 7I . 7 » 7 (m 2 r , — Q 2 , 0) but 
J r 7r * 7 * 7 (q|, q 2 , 0) is the relevant object to be used, where q 2 is to be integrated 
over. However, for large q\ the pion must be far off-shell, in which case the 
pion exchange effective representation becomes obsolete. Melnikov and Vain- 
shtein reanalyzed the problem by performing an operator product expansion 
(OPE) for q\~ q 2 ^> (<Zi + <Z2) 2 ~ rn 2 . In the chiral limit this analysis reveals 
that the external vertex is determined by the exactly known ABJ anomaly 
J r 7r77 (m 2 r , 0, 0) = — l/(47r 2 / w ). This means that in the chiral limit there is 
no VMD like damping at high energies at the external vertex. However, the 
absence of a damping in the chiral limit does not prove that there is no damp- 
ing in the real world with non-vanishing quark masses. In fact, the quark 
triangle- loop in this case provides a representation of the 7r°*7*7* amplitude 
given by 

F^, T {q 2 ,plpl) = (-4n 2 U)^. 7 *(q 2 ,plp 2 2 ) = 2m 2 C (m q ;q 2 ,p 2 ,p 2 ) 

2m 2 

[da] —= 2 2 2 ' ( 51 ) 

frig - a 2 a 3 p1 - a?,aiPi - aia 2 q z 

where [da] = daida2da3 (5(1 — a\ — a 2 — 03) and m q is a constituent quark 
mass (q = u, d, s). For p\ = p\ = q 2 = we obtain i^ » 7 « * (0, 0, 0) = 1, which 
is the proper ABJ anomaly. Note the symmetry of Co under permutations of 
the arguments (p 2 ,p 2 , q 2 )- For large p\ at p 2 ~ 0, q 2 ~ or p\ ~ p\ at q 2 ~ 
the asymptotic behavior is given by 

^ C o Q 7 M 7 (0,P?,0)^rln 2 r, (0,p 2 , Pl ) ~ 2 r Inr (52) 

m 2 

where r = The same behavior follows for q 2 ~ p\ at p\ ~ 0. Note that 
in all cases we have the same power behavior ~ ni 2 /p 2 modulo logarithms. 
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Table 3. LbL: Summary of most recent results for a M x 10 





no FF 


BPP 


HKS 


KN 


MV 


FJ 


TV°, 11,7]' 


+ O0 


85 ±13 


82.7 ±6.4 


83 ± 12 


114 ± 10 


88 ± 12 


axial vector 




2.5 ± 1.0 


1.7 ±0.0 




22 ± 5 


10 ± 4 


scalar 




-6.8 ±2.0 








0± 7 


TT, K loops 


-49.8 


-19 ±13 


-4.5 ±8.1 




0± 10 


-19 ±13 


quark loops 


62(3) 


21 ± 3 


9.7 ± 11.1 






21 ± 3 


total 




83 ±32 


89.6 ± 15.4 


80 ±40 


136 ± 25 


100 ± 39 



Thus at high energies the anomaly gets screened by chiral symmetry breaking 
effects. 

We therefore advocate to use consistently dressed form factors as inferred 
from the resonance Lagrangian approach. However, other effects which were 
first considered in [105] must be taken into account: 

1) the constraint on the twist four (l/g 4 )-term in the OPE requires hi = —10 
GeV 2 in the Knecht-Nyffeler from factor (|39"]l: Sa^ ~ ±5 ± 

2) the contributions from the /i and /{ isoscalar axial-vector mesons: Sa^ ~ 
±10 ± 4 (using dressed photons) 

3) for the remaining effects: scalars (Jo) ± dressed ix^ ,K^- loops ± dressed 
quark loops: Sa^ ~ ±2 ± 6 

Note that the remaining terms have been evaluated in [961 El! only. The split- 
ting into the different terms is model dependent and only the sum should be 
considered: the results read —5 ± 13 (BPP) and 5.2 ± 13.7 (HKS) and hence 
the true contribution remains unclean 11 !. 

An overview of results is presented in Table [3l The last column gives my 
estimates base on (96j [97l 11011 1105| . The "no FF" column shows results for 
undressed photons (no form factor). The constant WZW form factor yields 
a divergent result, applying a cut-off A one obtains [102] (a/7r) 3 C In 2 A, with 
an universal coefficient C = TV 2 ™ 2 /(48-7T 2 / 2 ); in the VMD dressed cases My 
represents the cut-off A — » My if My — > oo. 



5 Theory Confronting the Experiment 

The following Tab.S]collects the typical contributions to a p evaluated in terms 
of a determined via a e (|14j) . The world average experimental muon magnetic 
anomaly, dominated by the very precise BNL result, now is [11] 

o° x p = 1.16592080(63) x 10~ 3 (53) 



The problem seems to be the sizable negative scalar contribution of [97], which 
in |90| was estimated to be much smaller. Also the sign seems to be in question. 
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Table 4. The various types of contributions to a M in units 10 6 , ordered according 
to their size (L.O. lowest order, H.O. higher order, LbL. light-by-light) 



L.O. universal 


1161.409 73 


(0) 


e-loops 


6.194 57 


(0) 


H.O. universal 


-1.757 55 


(0) 


L.O. hadronic 


0.069 21 


(56) 


L.O. weak 


0.001 95 


(0) 


H.O. hadronic 


-0.001 00 


(2) 


LbL. hadronic 


0.001 00 


(39) 


r-loops 


0.000 43 


(0) 


H.O. weak 


-0.000 41 


(2) 


e+r-loops 


0.000 01 


(0) 


theory 


1165.917 93 


(68) 


experiment 


1165.920 80 


(63) 



World Ave 
Theory (g - * - 
Theory (t) 



3.2 cr 

1.3 cr 



— ISO — 50 O 



EW 1— loop 
EW 2— loop 
L.O. had 
H.O. had 
LbL 1995 
LbL 2001 
had T-r^^e 



lOO 200 
QED 



in units 
10 -io 



a u x lO 1 " — 11659000 



(relative uncertainty 5.4 x 10 7 ), which confronts the SM prediction 

af c = 1.16591793(68) x 1CT 3 . (54) 

Fig. [14] illustrates the improvement achieved by the BNL experiment. The 
theoretical predictions mainly differ by the L.O. hadronic effects, which also 
dominate the theoretical error. A deviation between theory and experiment 
of about 3 a was persisting since the first precise BNL result was released in 
2000, in spite of progress in theory and experiment since. 

Note that the experimental uncertainty is still statistics dominated. Thus 
just running the BNL experiment longer could have substantially improved 
the result. Originally the E821 goal was 8a e * p ~ 40 x 10~ n . Fig. [L31 illustrates 
the sensitivity to various contributions and how it developed in time. The 
dramatic (m M /m e ) 2 enhancement in the sensitivity of a^, relative to a e , to 
physics at scales M larger than m M , which is scaling like (m M /M) 2 , and the 
more than one order of magnitude improvement of the experimental accuracy 
has brought many SM effects into the focus of the interest. Not only are we 
testing now the 4-loop QED contribution, higher order hadronic VP effects, 
the infamous hadronic LbL contribution and the weak loops, we are reaching 
or limiting possible New Physics at a level of sensitivity which causes a lot of 
excitement. "New Physics" is displayed in the figure as the ppm deviation of 

«Ja M = a° xp - a^ c = (287 ± 91) x 10~ n (55) 

which is 3.2 a. We note that the theory error is somewhat larger than the 
experimental one. It is fully dominated by the uncertainty of the hadronic low 
energy cross-section data, which determine the hadronic vacuum polarization 
and, partially, by the uncertainty of the hadronic light-by-light scattering 
contribution. 
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CERN (79) 
KNO (85) Theory 
E821 (00) n+ 
E821 (01) n+ 
E821 (02) n+ 
E821 (04) it 
Average 
E969 goal 




208.0 ±6.3 



a^xio^-nesgooo 
I,, 



EJ 95 (e+e~) 

DEHZ03 | ^ V 

\ (+t) 
GJ03 (e+e~) 

SN03 (e+e- TH) i 
HMNT03 (e+e- inch) 
• (e+e') 

(+r) 
DEHZ06 (e+e") 
HMNT06 (e+e") 

ILbLBPP,HK,KN 
LbL FJ 
LbL M v 



TY04 



181.3 ± 16. 


[1.6 


a 


180.9 ±8.0 


[2.7 


a 


195.6 ±6.8 


[1.3 


a 


179.4 ±9.3 


[2.5 


a 


169.2 ±6.4 


[4.3 


a 


183.5 ±6.7 


[2.7 


a 


180.6 ±5.9 


[3.2 


a 


188.9 ±5.9 


2.2 


a 


180.5 ±5.6 


[3.3 


a 


180.4 ± 5.1 


[3.4 


a 


177.6 ±6.4 


[3.3 


n 


179.3 ±6.8 


[3.2 


a 


182.9 ±6.1 


2.9 


a 



Fig. 14. Comparison between theory and experiment. Results differ by different 
L.O. hadronic vacuum polarizations and variants of the LbL contribution. Some 
estimates include isospin rotated r-data (+r)). The last entry FJ06 also illustrates 
the effect of using different LbL estimations: 1) Bijnens, Pallante, Prades (BPP) [97] . 
Hayakawa, Kinoshita (HK) [98] and Knecht, Nyffeler (KN) [lOlj ; 2) my estimation 
based on the other evaluations; 3) the Melnikov, Vainshtein (MV) |105] estimate of 
the LbL contribution. EJ95 vs. FJ06 illustrates the improvement of the e+e~-data 
between 1995 and now (see also Tab.(2]|. E969 is a possible follow-up experiment of 
E821 proposed recently |115j 



As we notice, the enhanced sensitivity to "heavy" physics is somehow good 
news and bad news at the same time: the sensitivity to "New Physics" we are 
always hunting for at the end is enhanced due to 




by the mentioned mass ratio square, but at the same time also scale dependent 
SM effects are dramatically enhanced, and the hadronic ones are not easy to 
estimate with the desired precision. 
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Fig. 15. Sensitivity of g-2 experiments to various contributions. The increase in 
precision with the BNL g — 2 experiment is shown as a gray vertical band. New 
Physics is illustrated by the deviation (ajf p — aj l hc )/a° xp 

6 Prospects 

The BNL muon g — 2 experiment has determined as given by (|53[) . reaching 
the impressive precision of 0.54 ppm, a 14-fold improvement over the CERN 
experiment from 1976. Herewith, a new quality has been achieved in test- 
ing the SM and in limiting physics beyond it. The main achievements and 
problems are 

• a substantial improvement in testing CPT for muons, 

• a first confirmation of the fairly small weak contribution at the 2 — 3 a 
level, 

• the hadronic vacuum polarization contribution, obtained via experimental 
e + e~ annihilation data, limits the theoretical precision at the 1 a level, 

• now and for the future the hadronic light-by-light scattering contribution, 
which amounts to about 2er, is not far from being as important as the weak 
contribution; present calculations are model-dependent, and may become 
the limiting factor for future progress. 

At present a 3.2a deviation between theory and experiment is observed^ 
and the "missing piece" (|55|) could hint to new physics, but at the same time 
rules out big effects predicted by many possible extensions of the SM. 



It is the largest established deviation between theory and experiment in elec- 
troweak precision physics at present. 
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Fig. 16. Physics beyond the SM: leading SUSY contributions to g — 2 in supersym- 
metric extension of the SM 



Usually, new physics (NP) contributions are expected to produce contri- 
butions proportional to m^/M^p and thus are expected to be suppressed by 
Myy/M^p relative to the weak contribution. 

The most promising theoretical scenarios are supersymmetric (SUSY) ex- 
tensions of the SM, in particular the minimal one (MSSM). Each SM state X 
has an associated supersymmetric "sstate" X where sfermions are bosons and 
sbosons are fermions. This implements the fermion <-> boson supersymmetry. 
In addition, an anomaly free MSSM requires a second complex Higgs dou- 
blet, which means 4 additional scalars and their SUSY partners. Both Higgs 
holds exhibit a neutral scalar which aquire vacuum expectation values v\ and 
1)2- Typical supersymmetric contributions to a M stem from smuon-neutralino 
and sncutrino-chargino loops Fig. 1161 Some contributions are enhanced by 
tan/? = ^ which may be large (in some cases of order mt/mb ~ 40). One 
obtains [110j (for the extension to 2-loops see [111] ) 



susy • / sOt(M z ) (5 + tan 2 O w ) m n n ( „ 4a, m 

SUSY „ ( ) J V . 2 ^ -4 tan 1 In — 

487T sin W m A \ ir m M 



m = mgusY a typical SUSY loop mass and [i is the Higgsino mass term. In 
the large tan j3 regime we have 

|af SY | ^123x 10- 11 f 10 ° GeV Vtan/3. (57) 
1 ^ 1 \ m J 

a^ USY generally has the same sign as the /x-parameter. The deviation ([5 
requires positive sign(/i) and if identihed as a SUSY contribution 



m ~ (65.5 GeV) ^tan/? . (58) 

Negative /j, models give the opposite sign contribution to a M and are strongly 
disfavored. For tan /3 in the range 2 4- 40 one obtains 

m ~ 92 - 414 GeV , (59) 



precisely the range where SUSY particles are often expected. For a variety of 
non-SUSY extensions of the SM typically |a M (NP)| ~ C m^/M 2 where C = 
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0(1) [or 0{a/ir) if radiatively induced]. The current constraint suggests (very 
roughly) M ~ 1.7- 2.4 TeV [M ~ 87 - 121 GeV]. The C = 0(1) assumption 
is problematic, however, since no tree level contribution can be tolerated. 
For a more elaborate discussion and further references I refer to [112j . Note 
that the most natural leading contributions in extensions of the SM are 1- 
loop contributions similar to the leading weak effects or the leading MSSM 
contributions. However, mass limits set by LEP and Tevatron make it highly 
non-trivial to reconcile the observed deviation to many of the new physics 
scenarios. Only the tan (3 enhanced contributions in SUSY extensions of the 
SM for [i > and large enough tan j3 may explain the "missing contribution" . 
Two Higgs doublet models [113] have similar possibilities. Physics beyond the 
SM of course not only contributes to a M but also to other observables like 
to the branching fraction Br(b -> sj) = (3.40 ± 0.28) x 10~ 4 or to the W 
mass prediction M w = 80.392(29) GeV. In the imparity conserving MSSM 
the lightest neutralino is stable and therefore is a candidate for cold dark 
matter in the universe. From the precision mapping of the anisotropics in the 
cosmic microwave background, the WMAP collaboration has determined the 
relict density of cold dark matter to Slh 2 = 0.1126 ± 0.0081. This sets severe 
constraints on the SUSY parameter space (see for example [1 14J ) . 

Of course, for a specific model, one must check that the sign of the induced 
a^ p is in accord with experiment (i.e. it should be positive). 

Plans for a new g — 2 experiment exist |115j . In fact, the impressive 0.54 
ppm precision measurement by the E821collaboration at Brookhaven was still 
limited by statistical errors rather than by systematic ones. Therefore an 
upgrade of the experiment at Brookhaven or J-PARC (Japan) is supposed to 
be able to reach a precision of 0.2 ppm (Brookhaven) or 0.1 ppm (J-PARC). 

For the theory this poses a new challenge. It is clear that on the theory 
side, a reduction of the leading hadronic uncertainty is required, which actu- 
ally represents a big experimental challenge: one has to attempt cross-section 
measurements at the 1% level up to J/tp[T] energies (5[10] GeV). Such mea- 
surements would be crucial for the muon g — 2 as well as for a more precise 
determination of the running fine structure constant ckqed (E) . In particular, 
e + e~ low energy cross section measurements in the region between 1 and 2.5 
GeV |116l 1117] are able to substantially improve the accuracy of a)j ad ' 1 ' and 
a QED (M z ) [75]. 

New ideas are required to get less model-dependent estimations of the 
hadronic LbL contribution. Here, new high statistics experiments attempting 
to measure the 7r°7*7* form factor T{m^, —Q\, — Q|) f° r Q\ ~ Q\ an d a scan 
of the light-by-light off-shell amplitude via e + e~ — > e + e~7*7* — > e + e~77 
would be of great help. Certainly lattice QCD studies |118j will be able to 
shed light on these non-perturbative problems in future. 

In any case the muon g — 2 story is a beautiful example which illustrates 
the experience that the closer we look the more there is to see, but also the 
more difficult it gets to predict and interprete what we see. Even facing prob- 
lems to pin down precisely the hadronic effects, the achievements in the muon 
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g — 2 is a big triumph of science. Here all kinds of physics meet in one single 
number which is the result of a truly ingenious experiment. Only getting all 
details in all aspects correct makes this number a key quantity for testing 
our present theoretical framework in full depth. It is the result of tremendous 
efforts in theory and experiment and on the theory side has contributed a lot 
to push the development of new methods and tools such as computer alge- 
bra as well as high precision numerical methods which are indispensable to 
handle the complexity of hundreds to thousands of high dimensional integrals 
over singular integrands suffering from huge cancellations of huge numbers of 
terms. Astonishing that all this really works! 

Note added: After completion of this work a longer review article ap- 
peared [119] , which especially reviews the experimental aspects in much more 
depth than the present essay. For a recent reanalysis of the light-by-light 
contribution we refer the reader to [120] . which presents the new estimate 
a J; bL = (110 ±40) x 1CT 11 . 
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